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ABSTRACT 



o 
o 

In this paper, we give a description of a recent quantum algorithm created by Aharonov, Jones, and Landau 
for approximating the values of the Jones polynomial at roots of unity of the form e 27ri / k _ This description is 
given with two objectives in mind. The first is to describe the algorithm in such a way as to make explicit the 



underlying and inherent control structure. The second is to make this algorithm accessible to a larger audience. 
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1. INTRODUCTION 



In this paper, we give a description of a recent quantum algorithm created by Aharonov, Jones, and Landau 1 for 
approximating the values of the Jones polynomial at roots of unity of the form e 27 ™/ fc , for k a positive integer. 
We do so with two objectives in mind. The first is to describe the algorithm in such a way as to make explicit the 
underlying and inherent control structure. The second is to make this algorithm accessible to a larger audience. 

To avoid cluttering this exposition, we focus solely on the version of this quantum algorithm based on the 
Markov trace closure of a braid. The alternative platt closure version is left as an exercise for the reader. 

Readers already familiar with knot theory may want to skip to section 7. 



2. KNOT THEORY IS ... 

In its most general form, knot theory is the study of the fundamental problem of placement: 

The Placement Problem. When are two placements of a space X in a space Y the same or different? 



Placement Problem: Knot Theory | 

f \ M 3 or ff 3 ] 




Figure 1. The Placement Problem. 

In its most renowned form, knot theory is the study of the placement of a 1-sphere* .S^or a disjoint union of 
1-spheres) in real 3-space M 3 (or the 3-sphere S* 3 ), called the ambient space. In this case, "placement" usually 
means a smooth (or piecewise linear) embedding, i.e., a smooth homeomorphism into the ambient space. Such 
a placement is called a knot if a single 1-sphere is embedded ( or a link, if a disjoint union of many 1-spheres 
is embedded.) 

Two knots (or links) are said to be the same, i.e., of the same knot type, if there exists an orientation 
preserving autohomeomorphism^ of the ambient space carrying one knot into the other. Otherwise, they are said 
to be different, i.e., of different knot type. Such knots are frequently represented by a knot diagram, i.e., a 
planar 4-valent graph with vertices appropriately labelled as undercrossings/overcrossings, as shown in figure 2. 



*By 1-sphere we mean a circle. 

T A provenly equivalent definition is that two knots are of the same knot type if and only if there exists an isotopy of 
the ambient space that carries one knot onto the other. 
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Figure 2. A knot diagram. 

The fundamental problem of knot theory can now be stated as: 

The Fundamental Problem of Knot Theory. When are two knots of the same or of different knot type ? 

A useful knot theoretic research tool is Reidemeister's theorem, which makes use of the Reidemeister moves 
as defined in figure 3.: 

Theorem 2.1 (Reidemeister). Two knot diagrams represent the same knot type if and only if it is possible 
to transform one into the other by applying a finite sequence of Reidemeister moves. 




These are local moves ! 



Figure 3. The Reidemeister moves. 

The standard approach to attacking the fundamental problem of knot theory is to create knot invariants for 
distinquishing knots. By a knot invariant / we mean a map from knots to a specified mathematical domain 
which maps knots (or links) of the same type to the same mathematical object. Thus, if an invariant is found 
to be different on two knots (or links), then the two knots (or links) cannot be of the same knot typel 

The Jones polynomial is one such invariant, and indeed a very significant one. The Jones polynomial maps 
each knot (or link) to a Laurent polynomial with integer coefficients, i.e., an invariant with domain the Laurent 
ring Z [i, i -1 ] . If two knots (or links) have different Jones polynomials then they must be different. A definition 
of this famous knot invariant is given in section 6 of this paper. 

But before we define the Jones polynomial we will need to look at a group invented by Emil Artin, i.e., the 
braid group, which is now beginning to have an impact on quantum computing. 



3. THE BRAID GROUP B N 



Definition 3.1. The n-stranded braid group B n is is the group generated by the symbols b\, b%, . . . , b n —i 
subject to the defining relations 

bib i+1 b. L = b l+1 bib i+ i for 1 < i < n 
bibj = bjbi for \i - j\ > 2 

The braid group B n is easily understood in terms of diagrammatics, as shown in Figures 4 to 13. 

As illustrated in Figure 4, an element of the braid group B 3 can be thought of as a hatbox with the three 
points at the top connected to the three points at the bottom by smooth nonintersecting curves, called strands. 
As illustrated in Figures 5 and 6, two such elements of B n are equal if it is possible to continuously move within 
the hatbox the strands of one braid into the strands of the other without cutting or breaking the strands and 
without letting the strands pass through one another. Figure 5 shows two equal braids, and Figure 6 shows two 
braids that cannot be equal. As illustrated in Figure 7, the enclosing hatbox is usually omitted, but understood 
to be there. 

The product of two elements of B n is defined, as illustrated in Figure 8, by simply stacking one hatbox on 
top of the other. Under this definition of multiplication, it can be shown that every braid has a multiplicative 
inverse. An example is given in Figure 9 of the inverse of a braid. The braid group actually has a finite set of 
generators, as shown in Figure 10. A complete set of defining relations among these generators are shown in 
Figures 11 and 12. 




Figure 4. An example of a three 
stranded braid in B%. 



Figure 5. An example of two equal 
braids. 



I Two Unequal Braids | 





Figure 6. An example of two 
non-equal braids. 




Figure 7. Shorthand notation for 
braids. 



| Product of Braids | 




Figure 8. The product of braids. 
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| Inverse of oi a Braid | 



To construct the inverse of a braid, take the minor 
image of each crossing, and then reverse the order 
of the crossings 



Figure 9. The inverse of a braid. 



| Generators of the Braid Group B„ 
The braid group B n is generated by 
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Figure 10. The generators of the 
braid group. 



Figure 11. A diagrammatic 
illustration of one of the first set of 
defining relations. 
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Figure 12. A diagrammatic 
illustration of one of the second set 
of defining relations. 
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Figure 13. The closure (or trace) of 
a braid. 



4. TRANSFORMING BRAIDS INTO KNOTS AND VICE VERSA: THE CLOSURE 

(OR TRACE) OF A BRAID 

Why are braids of importance in knot theory? 

The answer to this questions begins by observing that every braid [3 can be converted into a knot (or link) 
by forming the closure (a.k.a., trace) (3 tr as shown in Figure 13. Then, of course, there is the famous Markov 
theorem telling us when two braids produce the same knot (or link): 

Theorem 4.1 (Markov). Two braids Pi and Pi produce the same knot (or link) under braid closure if and 
only if there exists a finite sequence of Markov moves that transforms one braid into the other. 

The two Markov moves M\ and Mi are are shown in Figures 14 and 15. 
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Figure 14. The Markov 1 move. 
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The Markov 2 Move 
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Figure 15. The Markov 2 move. 



Most amazingly, the process of transforming a braid into a knot (or link) can be reversed, as stated by 
Alexander's theorem: 



Theorem 4.2 (Alexander). Every knot (or link) is the closure of a braid. 



Returning to our original highly algebraic definition of the braid group B n , we should mention that each braid 
[3 in B n can be expressed as a product of the generators bi, 62, . . . , b n -i and their inverses 67/ 1 , b% 1 , • • • , b^li- 
Thus, each braid (3 can be written in the form 

0-TT& e(<) -b e(1) b e{2) ■■■6 £(/) 

i=l 

where e (i) = ±1 for i = 1, 2, . . . , I. We will call such a product in the generators &i, 62, • • • , &n-i an d their 
inverses b± , , . . . , a word defining the braid /3. Two words define the same braid (3 if and only if it is 
possible transform one into the other by applying a finite sequence of Tieze transformations. For a definition of 
Tieze transformations, we refer the reader to Crowell and Fox. 5 

Finally, we define the writhe of a braid /?, written Writhe {(3) , as the sum of the exponents in any word 
defining the braid. 

5. THE TEMPERLEY-LIEB ALGEBRA TL N (D) 

Our next stepping stone to the definition of the Jones polynomial is an algebra, called the Temperley-Lieb 
algebra. 

Definition 5.1. Let d be an indeterminate complex number. Then for each positive integer n, the Temperley- 
Lieb algebra TL n (d) is defined as the algebra with identity 1 generated by the identity 1 and the symbols 

Ei,E 2} . . . , E n _\ 

subject to the defining relations 

' E i E j =E j E i for >2 

< E i E i ± 1 E i =E> for 1 < i < n 

Ef = dEi for l<i<n 

The Temperley-Lieb algebra TL n (d) is easily understood in terms of diagrammatics*, as illustrated in Figures 
4 through 13. 

As shown in Figure 13, an element of the Temperley-Lieb algebra TL 3 (d) can be thought of as the algebra 
consisting of all linear combinations of rectangles, each rectangle with 3 points at the top and 3 points at the 
bottom connected by smooth non-intersecting curves, called strands. As shown in Figures 17 and 18, two such 
elements of TL n (d) are equal if it is possible to continuously move the strands of one into the strands of the 
other without cutting or breaking the strands and without letting the strands pass through each other. Figure 
17 shows two equal elements, and Figure 18 shows two non-equal elements. As shown in Figure 19, the enclosing 
rectangle is usually omitted, but understood to be there. 

The product of two elements of TL n (d) is defined, as shown in Figure 20, i.e., by stacking one rectangle 
on top of the other. As illustrated in Figure 21, if a circle should happen to arise as a result of the product, 
then it is simply replaced by the indeterminate d times the resulting rectangle with the circle omitted. The 
generators of the Temperley-Lieb algebra are shown in Figure 22. We leave, as an amusing exercise for the 
reader, the task of translating the complete set of defining relations for the Temperley-Leib algebra given above 
into diagrammatics. 

^The diagrammatic representation of the Temperler-Leib algebra is due to Kauffman. 18,20,21 




An element X of ILJd) 



Figure 16. 




Figure 17. 
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Figure 18. 
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Figure 19. 
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Figure 22. The Temperley-Lieb 
algebra TL„ (d) is generated by all 
formal sums of products of the 
above generators. 




Figure 23. The diagramatic 
definition of the Markov trace 

Tr n : TL n (d) — ► C. 



We would be amiss if we did not mention that there is a map Tr n : TL n (d) — > C from the Temperley-Lieb 
algebra TL n (d) to the complex numbers C, called the Markov trace, satisfying the following three conditions: 



• Tr n (1) = 1 

• Tr» (17) = Tr n (YX) for all X and Y in TL n (d) 

• If X e TL n (d), then Tr n+1 (XE n ) = \Tr n (X) 



A diagrammatic definition of the Markov trace is shown in figure 23. 



We will later need the following theorem: 

Theorem 5.2. The above three conditions uniquely determine the Markov trace, i.e., any map TL n (d) 
satisfying the above three conditions must be the same as that defined by figure 23. 



6. THE DEFINITION OF THE JONES POLYNOMIAL 



Let d be an indeterminate complex number, and let A also be an indeterminant complex number such that 
d = —A 2 — A~ 2 . Let TL n (d) be the corresponding Temperley-Lieb algebra, and let B n denote the n-stranded 
braid group. Then the Jones representation 

Pa ■ TL n (d) — ► B n 

is the group representation defined by 

hi i — ► AE. + A- 1 ! 
b- 1 i ^ A^Ek+Al 

where b\, b%, . . . , b n — l denote the generators of the braid group B n , and where 1, E\, E2, ■ ■ ■ , -EW-i denote the 
generators of the Temperley-Lieb algebra TL n (d) . 



We leave for the reader's amusement the exercise of verifying that the images under pa of the generators 
b\, 62, ... , b n -i satisfy the defining relations of the braid group. 

We are now finally in a position to define the Jones polynomial. 

Let (3 be an element of the n-stranded braid group B n , and let j3 Tr denote the knot (or link) constructed 
from the closure of the braid j3. Then the Jones polynomial Vm-r (t) of the knot (or link) j3 Tr is the Laurent 
polynomial in the polynomial ring Z [iji -1 ] over the integers Z given by 

Vprr (A~ 4 ) = -A 2Writhe( -^d n ~ 1 Tr n (p A (/?)) , 

where t — A^ A 1 where Writhe (/3 Tr ) denotes the writhe^ of the braid /3, and where Tr n (pa (/?)) denotes the 
Markov trace of the value of the Jones representation pa on the braid (3. 



7. THE REPRESENTATION $ : TL N (D) — ► CU (H N K ) OF THE TEMPERLEY-LIEB 

ALGEBRA TL N (D) 



Our objective is to describe the polytime quantum algorithm in AJL 1 for approximating values of the Jones 
polynomial Vqt t (t) at the primitive A:-th roots of unity t — e 2lxl l k for positive integers k. To this end, we begin 
by constructing a representation of the Temperley-Lieb algebra TL n (d) which carries the image of the Jones 
representation pa ■ B n — > TL n (d) onto a group of unitary transformations. 



Let Gfc denote the graph of k — 1 vertices and k — 2 edges given in Figure 24. 



The graph 
12 3 k-1 



Figure 24. The grapn-fi'fc ufk — °1 Jertices and k — 2 edges. 



The adjacency matrix of the graph Gk is easily seen to be 



M k = 



( 1 

1 1 

10 



V 



\ 




1 

1 0/ 



Moreover, d = 2cos(7r/fc) can be shown to be an eigenvalue of corresponding to the eigenvector 

A = (Ai, A 2 , . . . , A fe _i) = ( sin(n£/k) ) 1 < e<k 



We can now construct our representation as follows: 



Please refer to section 2 for a definition of writhe. 



Let V n ,k denote the set of paths in the graph Gk of length n starting at the vertex 1, and let H n ,k denote 
the Hilbert space with orthonormal basis 

{ \p) ■■ V G V n , k } 
with basis elements labelled by the paths p in V n .k- 

Interpreting as "to the left" and 1 as "to the right," we identify each length n path p with a binary string 
of length n. 

For each path p, let p 4_1 J be the subpath corresponding to the first i — 1 bits of p, let pL»— *+iJ denote the 
subpath corresponding to bits % up to and including bit z + 1 of p, and finally let p^ +2 denote the subpath of 
p corresponding to bits i + 2 up to and including the last n-th bit. Let (p) be the endpoint of the subpath 
p 1 -^ . Hence, e 4 (p) e {1, 2, . . . , k - 1}. 



Select d = 2 cos (n/k). Since A is related to d via the formula d = —A 2 — A 2 , we choose among the first 
four possible choices ± e l7r /( 2k ") for A the value 



A = e 



-»w/(2fc) 



We are now ready to define a representation 

<5 : TL n (d) — > C(7 (W n>fc ) 

of the Temperley-Lieb algebra TL n (d) into the group ring CU (Tt n ,k) of the group U (H n ,k) of unitary transfor- 
mations on the Hilbert space H n ,k by specifying the images 

$i = *(£i) 

of each of the generators £?j of the Temperley-Lieb algebra TL n (d), taking great care to make sure that the $i's 
satisfy the same defining relations as the E^s. 

We define <j>j as: 



r 



®i \P) = < 



A ej(p)-1 _|_ \/ A ei(P>- lA °i(PH 



A »i(P> 



A =i(P> 



if pLi-i+lJ = 00 



1^-^10^+2} if pLi-i+lJ = 01 



I o 



|pi-lJ01pL<+2\ + ^i(s)±i|p) if pLi-i+iJ = io 

if pLi-i+lJ = 11 



We leave for the reader the exercise of showing that the transformations <frj satisfy the defining identities of the 
2?i's, i e-, the identities 

= for |i - j| > 2 

$i$i±i$i = * 4 for 1 < i < n 



for 1 < i < n 

and hence that $ is a legitimate representation of the Temperley-Leib algebra TL n (d). 



8. CONSTRUCTING A TRACE TR COMPATIBLE WITH THE MARKOV TRACE 

TR N 



We next need to construct a trace Tr on the image of the representation <I> : TL n (d) — ► CU (7~L n ,k) which is 
compatible with the Markov trace Tr n , i.e., a trace Tr such that the following diagram is commutative 

TL n (d) Im($)cCC/(W„, fe ) 
Tr n \ I Tr 

C 

For this construction, we need the following lemma: 

Lemma 8.1. The representation <£> : TL n (d) — ► CU (W n ,fe) maps each ket \p) to a linear combination of kets 
each labeled by a path of the same length as the path p, and each having the same endpoint as p. 

An immediate corollary is: 

Corollary 8.2. Let V n .k,m be the subset of V n ,k of all paths p in Gk of length n that start at the vertex 
1 and end at the vertex m, where 1 < m < k, and let H n ,k,m be the Hilbert subspace of H n ,k defined by the 
orthonormal basis {\p} : p G V n ,k,m}- Then the representation $ : TL n (d) — ► CU (H n ,k) splits into the direct 
sum of representations 

fc-i 

$ = 0$(™) , 

m— 1 

fc-1 

where $( m ) : TL n (d) — > CU (H n ,k,m) is the representation arising from the projection ^^CU (H n ,k,m) — f 

m—l 

CU (H n ,k,m)- Hence, the image Im(3>) of the representation $ : TL n (d) — > CU (H n ,k) lies in the direct sum 
of the algebras CU {H n ,k,m), 1 < m < k, i.e., 

fe-i 



Im($) C 0C£/(W n , fc , m ) 



The above corollary gives us the latitude of searching for a compatible trace from among all the traces 
Tr : Im ($) — ► C which are constructed by taking any linear combination of the standard traces Tr m : 
CU (H n ,k,m) — * C Of these traces, the desired compatible trace is found to be the one given in the following 
theorem: 

THEOREM 8.3. Let \i = sin^/fc) be the components of the eigenvector A given in section 6, and let N = 
J2 e Z 1 Xidim(H n ,k,e)- Then the trace Tr : Im($) — > C defined by 

^ fe-i 

Tr =N^ ¥n ' 
i=i 

is compatible with the Markov trace, i.e., Tr is a trace such that the diagram 

fe-i 

TL n (d) Im ($) C 0CC/ (H n , k ,e) C CU (H n , k ) 

Tr n \ | Tr 

C 



is commutative. In other words, Tr is a trace such that Tr n = Tr o $. 



Since the Markov trace Tr n : TL n (d) — ► C is the unique trace satisfying the following three conditions 

• Tr n (1) = 1 

• Tr n (XY) = Tr n (YX) for all X and Y in TL n (d) 

• If X e TL n (d), then Tr n+1 (XE n ) = \Tr n (X) 

the proof of the above theorem consists simply in verifying that Tr o $ : TL„ (d) — > C satisfies each of these 
conditions. 



9. INTERMEDIATE SUMMARY 

But where are we in regard to our objective of creating a quantum algorithm for approximating the value of the 
Jones polynomial at a root of unity of the form e 27 "/ fc , where k is an arbitrary positive integer? 

For a knot (or link) given by the closure [3 Tr of an n-stranded braid (3, we have seen that the Jones polynomial 
is given by the expression 

V P tt (t) = -A 2Wnth <^d n - x Tr n (p A (/?)) , 
where d and A are indeterminate complex numbers related by the equation d = —A 2 — A~ 2 , and where t = A~ A . 

Setting A = e - 2 W 2 fc ( wn ich implies d — 2 cos (n/k) and t = e 2m / k ), we have the value of the Jones polynomial 
at t = e 2 W fc i s given by 

Vprr (e 2 "/ fe ) = -A 2Wnthe< >P ) d n - 1 Tr n (p A (/?)) . 

Since — A 2Wrlthe ^ d n ~ x is easily computed, the task of determining Vqtt ( e 27r4 / fc ^ reduces to that of evaluating 
the trace 

Tr n { PA {(3)) . 

But from the previous two sections, we have 
Tr n ( PA (/?)) = Tr n [(* o p A ) (/?)] = Tr- 

where Tr denotes the standard trace, where A m = sin (irm/k), and where N — Y^m=\ dim (H n ,k,m)- Thus, 
our objective reduces to finding the trace of each of the following k — 1 unitary transformations (called global 
gates) 

U(m) = ^p(m) o (/3) , 1 < m < k 

If the knot (or link) is given by the closure of a braid (3 defined by the word 

6-T\b eW -fo e(1 V (2) ...& e(i) 




fc-i 



(m) 



° PA (/?) 



-^A ro Tr[($No M ) 09) 



where b\, 62, ... , b n -\ are the generators of the braid group B n , and where e (i) — ±1 for i — 1, 2, . . . , L, then 
each unitary transformation — (<£>( m ) o p^j (/?) is given by 



U im) = Y[ (u 



Am) 



where [/■ m - ) denotes the unitary transformation (called an intermediate gate) 



U. 



(m) 



($ W oPa) (&*) , l<m<fc, 1<j < 



Thus, the trace Tr n (p^ (/?)) we seek to approximate is given by the following expression 



fe-i 



Tr„ ( PA (/?)) = - £ A m 7Y C/(™) = - £ A m Tr 



AT 



n K?>>) 



We are finally in a position to describe the quantum algorithm found in 1 for approximating the Jones polyno- 
mial VpTT (t) at t = e 2jri / fc as a quantum algorithm consisting of the completion of two sequences of steps, called 
phases. The first is a preliminary phase called the compilation phase. After completion, the compilation 
phase is immediately followed by a second phase, called the execution phase. 



10. THE COMPILATION PHASE 



The compilation phase is illustrated in figure 25. 
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Figure 25. The compilation phase. 



It consists of the following steps: 



Software Compilation. On receiving a regular diagram of a knot (or link) K as input, the algorithm described 
by Alexander's theorem is executed to produce a regular diagram of a braid [3 whose closure gives the knot (or 
link) K. (See Birman. 4 ) An algorithm, called braid combing, is then applied to the planar diagram of the 
braid /?, producing as output a word 

1=1 

describing the braid, and also as a side effect, producing as output the integer n giving the number of strands in 
(3. (Once again, see Birman. 4 ) This word can be thought of as a computer program which will later be compiled 
into hardware. 



First Hardware Compilation. Upon receiving as input the integers k and n, use the Kitaev-Solovay, 281 
theorem to implement (translate into hardware) good approximations of each intermediate gate [7 ■ , 1 < j < n, 
1 < m < k as a product of polynomially many elementary gates. 



Second Hardware Compilation. Use the braid word 



= 1 



to implement (i.e., to physically construct) the global gates [/( m ) 1 < m < k from the intermediate gates U, 
In other words, construct [/( m ) using the formula 

<"-•= n ( c «)"" ■ 



(m) 



Third Hardware Compilation. For each m (1 < m < k), construct from the global gate [/( m ) two quantum 
subroutines QRe m and QIm m , which upon input \p) (p <G V n ,k,m), output at random a '0' or a '1' according to 
the following probability distributions 

with probability Prob (0) = \ + \ Re (p\ \p) 

1 with probability Prob (1) = \ - \ Re (p\U^ \p) 

with probability Prob (0) = \ - \ Im (p\ U {m ^ \p) 

1 with probability Prob (I) = \ + \ Im (p\U (m) \p) 

where Re (p\U^ \p) and Im (Tp\U^ m ^ \p) denote respectively the real and the imaginary parts of the bracket 
{p\U^\p). Thus, if QRe m (\p)) and QIm m (\p)) arc repeated executed, then we obtain respectively approxi- 
mations of the real and imaginary parts of the bracket {p\U^ m ^\p), as displayed below: 

(#0' S -#l' S )/(#0',s + #l' S )«Re<p|C/( m )|p) for QRe m (|p» 
(#l'.s - #0'«) / (#0's + #1'*) w lm(p\U^\p) for QIm m (\p)) 



QRe m (\p)) 



QIm m (\p)) 



These two quantum subroutines QRe m and QIm m can be implemented using what has now come to be 
known as the Hadamard test. A wiring diagram defining the quantum subroutine QRe m (\p)) is given in 
figure 26. Input to QRe m consists of the state \p) and an ancillary qubit in state |0). After execution of the the 
wiring diagram, measurement of the ancillary qubit with respect to the standard basis will produce the desired 
probability distribution. Similarly, Figure 27 defines the quantum subroutine QIm m . 



Computing the Trace with Quantum Mechanics 
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Figure 26. A wiring diagram describing the quantum subroutine QRe m . Input consists of the 

state \p) and an ancillarary qubit in state |0). The Hadamard gate is denoted by H. After 
execution of the first hadard gate, the ancillary qubit is used to control the global gate U^ m \ 
Measurement of the ancillary qubit after the execution of the second Hadamard gate produces 

the desired probability distribution. 



Computing the Trace with Quantum Mechanics 



|o). 



\±cb(0)=±~Im(p\U\p) 



y 



J*ob(l)=± + hm{p\U\p} 

Figure 27. A wiring diagram describing the quantum subroutine QIm m . This subroutine is the 
same as the quantum subroutine QRe m except for the additional single qubit phase gate 
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11. THE EXECUTION PHASE 



The compilation phase is then followed by the execution phase as described by the pseudocode given below: 



Execution Phase for Estimating V^rr 



Algorithm A JK (n, k) 
Trace = 

LOOP m = 1 . . . k - 1 

Trace m — AND A m = sin (irm/k) 
LOOP p € "P„, fc , m 

i?e = AND 7m = 

LOOP Iteration = 1 . . . Numberof Iterations 

RealBit = QRe. m (p) and ImgBit = QIm m (p) 

Re = Re + {-l) Realmt and Im = Im - (-i) Im ^t 

END Iteration-LOOP 

Trace m = Trace m + (Re + i * Im) / Numberof Iterations 

END p-LOOP 

Trace = Trace + X m * Trace m 
END m-LOOP 
output Trace 
end Algorithm AJK 



This phase consists of three nested loops. The innermost loop calls the quantum subroutines QRe m and 
QIm m . The parameter Numberof Iterations is chosen according to the Chernoof-Hocffding bound to provide 
the desired accuracy for the approximation. If resources are available, the outermost iteration loop can be 
replaced by a parallel implementation. 



12. CONCLUSION 

Indeed, much more could be said about the AJL quantum algorithm for the Jones polynomial. But, for the 
time being, we will leave that task to one of our future forthcoming papers on this subject. 
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